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STREAMWISE VORTICITY GENERATION IN LAMINAR AND TURBULENT JETS 

AYODEJI O. DEMUREN* AND ROBERT V. WILSONt 

Abstract. Complex streamwise vorticity fields are observed in the evolution of non-circular jets. Gen- 
eration mechanisms are investigated via Reynolds-averaged (RANS), large-eddy (LES) and direct numerical 
(DNS) simulations of laminar and turbulent rectangular jets. Complex vortex interactions are found in 
DNS of laminar jets, but axis-switching is observed only when a single instability mode is present in the 
incoming mixing layer. With several modes present, the structures are not coherent and no axis-switching 
occurs. RANS computations also produce no axis-switching. On the other hand, LES of high Reynolds 
number turbulent jets produce axis-switching even for cases with several instability modes in the mixing 
layer. Analysis of the source terms of the mean streamwise vorticity equation through post-processing of the 
instantaneous results shows that complex interactions of gradients of the normal and shear Reynolds stresses 
are responsible for the generation of streamwise vorticity which leads to axis-switching. RANS computations 
confirm these results, k — ( turbulence model computations fail to reproduce the phenomenon, whereas 
algebraic Reynolds stress model (ASM) computations, in which the secondary normal and shear stresses are 
computed explicitly, succeeded in reproducing the phenomenon accurately. 

Key words, jets, DNS, LES, Navier-Stokes, turbulence modeling 

Subject classification. Fluid Mechanics 

1. Introduction. Experiments ([3, 4]) have shown that three-dimensional (3-D) jets can be used to 
enhance mixing and entrainment rates in comparison to axi-symmetric jets. A fundamental understanding 
of the dynamics of complex, turbulent jets is required for their prediction and control. 

Zaman [7] used streamwise and azimuthal vorticity dynamics to explain the presence or absence of axis- 
switching in experimental measurements of 3:1 aspect ratio rectangular jets with different, initial conditions. 
This study showed that the presence of streamwise vorticity pairs with outflow rotation (pumping fluid from 
the core to the ambient perpendicular to the major axis plane) produced axis switching while pairs with 
the opposite sense of rotation did not. However, in jets with no streamwise vorticity at discharge some 
other mechanism must originate it. Hussain and Husain [2] have explained that vortex self-induction which 
is an inviscid mechanism could be responsible. In the current study, numerical simulations of laminar and 
turbulent rectangular jets are performed to investigate the origin of streamwise vorticity and the mechanism 
for axis-switching through direct computations and analysis of terms in the streamwise vorticity equation. 

2. Mathematical Formulation. The partial differential equations governing the incompressible jet 
fluid flow are variants of the Navier-Stokes equations which can be written in Cartesian tensor form, for 
dimensionless variables as: 


(2.1) gg i+u ®h. = dP | 1 ^Uj dr,j 

dt 3 dxj dx, Re D dxjdx.j dxj 
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and 


where u, represents components of the velocity which are resolved in the computations, and r,j = uJuJ-UiUj 
is the Reynolds stress which must be modeled in terms of the resolved velocity field. At low Reynolds 
numbers, all scales of the velocity are resolved, therefore no modeling is required and the Reynolds stress is 

zero. Two variants are considered for high Reynolds number jet flows. 

In the LES formulation, the larger scales of the flow are resolved, and the smaller scales are modeled 
in terms of sub-grid scale (SGS) Reynolds stress. The Smagorinsky eddy-viscosity model is utilized in the 
present study to approximate the SGS Reynolds stress as: 


(2.3) 


nj - 6 fr kk = -2C s A 2 |S|Sy 


where <5^ is the Kronecker delta, A is a (dimensionless) leng thjicale a ssociated with the grid size, 5<j = 
+ is the resolved strain rate tensor, and \S\ = y/2 . The model coefficient, C s , is set 

to the constant value of 0.01. The budget of the time-averaged momentum equations show ([5]) that, in 
present model studies, magnitudes of the SGS stresses are only a small fraction of magnitudes of the resolved 
stresses, hence the use of a more sophisticated SGS model is unwarranted. 

In the RANS formulation, only the mean flow is resolved. Both small and large scales are time-averaged, 
and their effects are modeled via the turbulent Reynolds stress. Two turbulence models are considered 
for approximating the Reynolds stress; the k - e model which is based on a Boussinesq eddy viscosity 
hypothesis, similar to Eq. 2.3, and an algebraic Reynolds stress model (ASM), which is the simplest form 
for a second-moment turbulence closure. In the former, the Reynolds stress is calculated from: 


(2.4) T a - ~$ Tkk - 

In the latter, it can be shown ([1]) that the Reynolds stress has two parts; a linear part with exactly the 
same form as Eq. 2.4, and a non-linear part. Thus, 


(2.5) Aj = r * + T lj 

where r* is given by Eq. 2.4 and, for the secondary Reynolds stress components, which are critical for the 
generation of streamwise vorticity, j has the form: 


( 2 . 6 ) 


r'ij = k 


2«i % - a ?T 


dUi dUi 

dxi dxj 


In Eq. 2.6, T = k/e is the turbulent time scale, and a\ and e *2 are model coefficients. Hence, the 
second term of this equation, which is non-linear in T also represents a ratio of turbulent to mean time 
scales. In addition, it embodies the anisotropy of the turbulent field. At low Reynolds numbers, the ratio 
will be small, and the contribution of the term will tend to zero. This ASM, though highly simplified was 
shown by Demuren and Rodi [1] to be capable of successfully predicting turbulence-driven secondary flow in 

non-circular ducts. 
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For DNS or LES, Eq. 2.1 is .discretized temporally with explicit Runge-Kutta (RK) schemes and spatially 
with implicit compact finite difference schemes. Detailed descriptions are given in Wilson [5], For R ANS. 
the commercial code CFX-TASCflow, which uses second-order spatial discretization, is utilized. 


3. Model Problems. Spatial simulations of jet flows are performed in this study in which a fixed region 
of the flow is computed and disturbances grow' in the streamwise direction. Rectangular jets are simulated 
wdiich have a nominal aspect ratio of 2:1. Reynolds numbers based on the core velocity and the equivalent, 
diameter (D) are 750 for laminar cases and 75,000 for turbulent cases. Details of the DNS and LES can be 
found in Wilson and Demuren [6]. Typical computations utilized (129 x 129 x 129) grid nodes for a domain 
size of (12 x 10 x 10) and required several hours of CPU time on the CRAY C-90 supercomputer. The RANS 
simulations were performed for a one-quarter segment, taking advantage of symmetry conditions. For this 
quarter-segment, the domain size relative to the jet diameter was (12 x 3 x 3) and the grid distribution was 
(80 x 40 x 40). Convergence was obtained in about 200 iterations which required about 6 hours of CPU time 
on a SUN ULTR A 10/360 workstation. 


4. Results and Discussion. Figure 4.1 shows contours of the instantaneous total vorticity for rect- 
angular jets at low Reynolds number (=750) along major and minor axes, respectively. The jet in the top 
frames (a, b) has fundamental mode instability imposed in the mixing layer at its inlet. This would be typical 
of a jet issuing from a pipe with a laminar boundary layer. The lower jet (c, d) has broad mode instabilities 
imposed in the mixing layer at its inlet, which would be typical of a jet issuing from a pipe with a turbulent 
boundary layer. The former shows well-organized structures, which by the end of the simulation at 10 diam- 
eters had led to a shrinking of the jet in the major-axis plane and an expansion in the minor axis plane. In 
the latter, the structures are not so well-organized and no axis-switching can be discerned. This result would 
also be typical of a natural unforced jet. Time-averaged results of the streamwise vorticity and velocity, 
shown in Figure 4.2 confirm that, in the mean, no axis-switching occurs. The structures are organized, in 
the mean, but they lead to a gradual evolution of the jet cross-section from rectangular at inlet to circular 
in the far-field. The streamwise vorticity field is such that would produce this gradual transition. Therefore, 
at low Reynolds number, a natural rectangular jet or one with a broad mode of instabilities in its mixing 
layer would not experience the phenomenon of axis-switching, whereas a rectangular jet with a fundamental 
instability mode forcing would. It was also shown in Wilson and Demuren [6] that with the addition of the 
sub-harmonic mode the tendency is towards jet bifurcation. Figure 4.3 shows the instantaneous total and 
streamwise vorticity from the LES of rectangular jets at high Reynolds number (=75,000), with broad-mode 
instabilities, along minor and major axes, respectively. There is no discernible streamwise vorticity in the 
first 2 diameters, and the evolution trends are quite different from those of the laminar jet: there is expansion 
in the minor-axis plane and contraction in the major-axis plane. Corresponding time-averaged results of the 
streamwise vorticity and velocity, shown in Figure 4.4. So what is the origin of axis-switching, and in what 
flow situations would it be expected to occur? 


4.1. Mechanism for Streamwise Vorticity Generation. The mechanism for streamwise vorticity 
generation can be examined by considering the time-averaged streamwise vorticity equation. Following 
Demuren and Rodi [1], the streamwise vorticity equation has the form: 
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Fig. 4.1. Instantaneous total vorticity for DNS of rectangular jet at Re = 750, with (a, b) fundamental and ( c , d) broad 
mode instabilities . 


(4.1) 


d<Ji dUi _ d / duiU2 _ 9 2 

Ul 9x7 ~ Ui 9x7 ~ 9x7 V 9x 3 9x 2 J 9x 2 9x 3 

^3 

_9^_ _9^_\ 1 f 9 V \ 

dx\ Ox 1 , ) 2 ' Re D \dxjdxj ) 

1 V V — V- 

v4s ^6 




The components of the vorticity vector are given by: 

dU 2 dU-i _ 9f7i 9t / 3 _ 9t / 2 _ 9t/i 

9x 3 9x 2 ’ 2 9x 3 dx\ ’ 9xi 9x 2 

The „4i terms represent the convection of streamwise vorticity by the mean velocity while the .4 2 terms 
represent the tilting and stretching of the vorticity vector by gradients of the mean velocity. Terms .4 3 , 
A 4 , and .4 r , contain the turbulent stresses and act to produce or destroy streamwise vorticity. In particular, 
terms A 4 and A 5 contain the effect of the turbulent normal and shear stresses, respectively. The diffusion 
of streamwise vorticity is given by the terms, A $ . 

The generation of streamwise vorticity through vortex stretching and tilting (terms A 2 ) is known as 
secondary motion of Prandtl’s first kind. This is an inviscid mechanism which may be present in laminar 
as well as turbulent flows. Generation of streamwise vorticity by gradients of the turbulent stresses (terms 
.4 3 , . 44 , and .4 5 ) is known as secondary motion of Prandtl’s second kind. Since there is no mean streamwise 
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Fig. 4.2. Cross-sectional contours of mean streamwise vorticity and velocity for DNS of rectangular jet at Re = 750 with 
broad mode instabilities. 


vorticity imposed at the inflow of the current simulations , its generation downstream of the inflow can only 
be explained by tilting and stretching of the azimuthal vorticity and/or by gradients of the turbulent normal 
and shear stresses. 

The terms of the mean streamwise vorticity equation were computed by post- processing the instantaneous 
results from the DNS and LES of the spatially-developing rectangular jet with broad mode forcing. An 
examination of these can illuminate the origin of streamwise vorticity in these jets. The results show that 
the streamwise vorticity component in the first two diameters is several orders of magnitude smaller than 
the azimuthal component at the inflow. Streamwise vorticity undergoes rapid growth in the region 2.5 < 
x/D < 4. 

In the presence of fundamental mode instability, rollers are formed in the mixing layer of the jet, which 
would be rectangular vortex rings in the case of a rectangular jet. But from the Biot-Savart law, vortex 
rings have induced velocities which are proportional to their curvature. Hence, in the corners or along 
the minor-axis sides, the induced velocities would be larger than those along the major-axis sides. The 
differential induced velocity lead to deformation of the vortex rings and hence the jet, shown in Figure 4.1. 
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Fig. 4.3. Instantaneous total and streamxvise vorticity for LES of rectangular jet at Re — 75,000 with broad mode 
instabilities a) minor axis and b) major axis. 


The mechanism is inviscid and the jet deformation would proceed in this way, whether the jet is laminar 
or turbulent, so long as the vortex rings represent the dominant structure in the jet. This has been widely 
observed in the literature, and indeed, in many experimental and numerical studies, very strong forcing at 
the fundamental mode is applied not only to the mixing layer, but also to the whole jet. The consequence 
is that the jet evolution is largely governed by this discrete-mode forcing. 

In natural jets or jets with broad-mode forcing, organized vortex rings, if present at all, play only a 
minor role. In the laminar jet, the A 2 term is dominant and shows the gradual transfer of vorticity from 
the azimuthal components to the streamwise component. This process continues to increase the streamwise 
component, even beyond xfD = 8.25. The jet shape distorts from its initial rectangular to a circular one 
in the process. No axis-switching occurs in this case. In the turbulent jet, on the other hand, streamwise 
vorticity growth is initially due to the A 2 and A 4 terms. The A 5 term quickly builds and surpasses the A 2 
term, which subsequently stagnates and decays. Both the A 4 and A 5 terms grow to reach their maximum 
value before the end of the potential core. 

To obtain a global picture of the interaction of the source terms of the streamwise vorticity equation, 
detailed cross-sectional plots are required (not shown here) to analyze the distribution of the source terms 
within the cross-sectional plane and to explain the generation of the streamwise vorticity and axis switching 
mechanism. Cross-sectional contours of the A 2 , A 4 , A 5 , and the A 4 + A 5 terms, at x/D = 3, along with 
the mean streamwise vorticity and velocity show that eight pairs of streamwise vortices are present and 
their locations and signs are consistent with that of the A 4 + A5 term. Hence, it is shown that gradients 
of the difference of the turbulent normal stresses are responsible for the initial generation of streamwise 
vorticity. Contour patterns of streamwise vorticity and of the A 4 term show little bias towards the major 
or minor axis planes, while contours of the A 4 4- A5 term begin to show an anisotropic pattern. Further 
downstream at x/D = 3.75, contours of streamwise vorticity reveal that the vortex pairs orientated about 
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Fig. 4.4. Cross-sectional contours of mean streamwise vorticity and velocity for LES of rectangular jet at Re = 75,000. 


the minor axis plane are of greater strength than those along the major axis plane. Contours of the term 
show a coalescence of the like-signed contours present at x/D = 3 along the major axis plane, resulting in 
a weakening of the opposite-signed patterns of the -44 term along the major axis plane. This is apparent 
from the contours of the ,44 4- ,4 5 terms and explains the directional bias present in the streamwise vorticity 
contours. The sense of rotation of the vortex pairs about the minor axis plane is such that core velocity fluid is 
pumped from the jet centerline to the ambient leading to a distortion of the initially rectangular mixing layer. 
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The four pair of streamwise vortices oriented along the minor axis strengthen until the end of the potential 
core where the initially flat major axis sides develop a bulge and the corners are flattened, leading to the 
observed diamond shape at x/D = 5.25 (Figure 4.4) and the eventual switching of the major and minor axes. 
Comparison of the structures of the A 4 , A 5 terms and the streamwise vorticity suggests that the -4 4 term 
acts largely as a source and the A$ term as a sink, and the difference between them produces the streamwise 
vorticity. The inviscid mechanism, based on the Biot-Savart. law or .4 2 plays no significant role. These are 
exactly the roles that Demuren and Rodi [1] found that the turbulent terms played in generating streamwise 
vorticity in non-circular ducts. This similarity suggests that RANS computations based on second-moment 
closure models should reproduce these effects whereas those based on Boussinesq-type eddy-viscosity models 
should fail. 



Fig. 4.5. Jet half-widths for RANS calcAtlation of rectangular jet at Re = 75,000; ASM solid lines , k - c dashed lines . 

Results of RANS computations of the turbulent rectangular jet at a Reynolds number of 75,000 are 
presented in Figures 4.5 and 4.7. Figure 4.5 shows the growth of the jet half-widths. The ASM results 
exhibit growth along the minor axis and decline along the major axis, with a cross-over between the major 
and minor axes at x/D of 3. On the other hand, the k - e model results show growth along both axes, but 
at a faster rate along the minor axis so that by the end of the calculation domain at x/D of 12, both widths 
tire nearly equal. Flow patterns are shown in Figure 4.6 for three cross-stream locations, x/D of 1.25, 2.5 
and 5. The k-t model results show entrainment of ambient fluid from both sides and some transfer of the 
fluid from the narrow side to the wide side which leads to faster growth of the jet along the minor axis, 
and the approach towards a circular cross-section. This result mirrors that shown in Figure 4.2 for laminar 
rectangular jets rather than the turbulent jet ones of Figure 4.4. The ASM results show much more complex 
flow patterns. At x/D of 1.25, the potential core is largely undisturbed, but the mixing layer region is highly 
distorted. There is a strong reverse flow region, and the sense of the circulation is to move fluid from the 
narrow side to the wide side which produces strong distortion of the streamwise velocity contours. By 2.5 
diameters, the jet shape is now highly distorted, though a potential core still exists which maintains its 
rectangular shape. But by 5 diameters, at the end of the potential core, the major and minor axes of the jet 
have switched, and there is now one major circulation which sense is to continue the axis-switching process. 
Figure 4.7 shows the streamwise vorticity contours. At x/D of 1.25, there are two pairs of vortices, consistent 
with the LES results of Figure 4.4. The vortex strength increases up to 2.5 diameters, where only one pair 
remains. The bulk of the jet distortion occurs in this region. By 5 diameters, the strength of the vortex 
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Fig. 4.6. Mean velocity vectors and streamxuise velocity contours at X/D = 1.25,2.5,5.0, for HANS computations of 
rectangular jet at Re = 75,000; ASM - left side , k - < - right side. 

has decreased by almost one order of magnitude below the peak value. What is seen are the final stages of 
the axis-switching process. The LES and RANS computations with ASM show good qualitative agreement, 
though the vortex strengths appear to be higher and the axes appeared to switch earlier in the latter. Jet 
evolution has been shown in the present study, and in many experiments, to be strongly influenced by the 
instabilities or turbulence in the inlet mixing layer, so quantitative comparisons would require a closer match 
of these conditions. Furthermore, the version of second-moment closure utilized here is highly simplified and 
was found to over-predict secondary flow generation in some non-circular ducts, though essential features 
are reproduced. It w*as adopted here for its simplicity, as the lowest level of modeling required to confirm 
the explanations for the origin of mean streamwise vorticity found in the analysis of the LES results. 

5. Concluding Remarks. Three-dimensional simulations of laminar and turbulent jets with rectangu- 
lar cross-section were performed. At low Reynolds numbers DNS were performed, while at higher Reynolds 
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Fig. 4.7. Mean velocity vectors and. streamwise vorticity contours at X/D — 1.25,2.5,5.0, for RANS computations of 
rectangular jet at Re = 75,000; ASM - left side , k - e - right side . 


numbers LES or RANS computations were performed. Origin of secondary flow or streamwise vorticity under 
different conditions of the spectral content of the initial jet mixing layer or Reynolds number are investigated 
through budgets of the mean streamwise vorticity derived from the DNS and LES data. Inviscid mecha- 
nisms, which can be explained by the evolution of vortex rings are responsible for the distortion of laminar 
jets with discrete mode forcing. Natural laminar jets evolve in an uneventful manner from rectangular to 
circular cross-section. In turbulent jets, streamwise vorticity is generated by terms involving derivatives of 
the secondary Reynolds normal and shear stresses. Inviscid mechanisms play no role. These results are 
confirmed by RANS computations with ASM. 
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